The absence of surface states for LiFeAs 
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We investigate the cleaving behavior of LiFeAs and determine its surface electronic structure by 
detailed density functional calculations. We show that due to the neutral surface of LiFeAs after 
cleaving, barely any influence of the surface on the electronic states is present. Therefore the data 
of surface sensitive probes such as angle resolved photoemission spectroscopy (ARPES) represent 
to a high degree the bulk electronic structure. This we highlighted by a direct comparison of the 
calculations to ARPES spectra. 
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The use of angle-resolved photoemission (ARPES) and 
scanning tunneling spectroscopies to study the electronic 
structure of iron pnictides can turn out to be crucial to 
understand the superconductivity in these materials. On 
one hand, both methods are very informative by provid- 
ing the detailed distribution of electron density in real 
and reciprocal spaces, but on the other hand it is not 
clear whether the obtained information is related to the 
bulk or surface. 

In a recent paper [l[ we found a pronounced sur- 
face electronic structure in LaOFcAs, which appears 
to be present in published angle-resolved photoemission 
(ARPES) data. 0-01 LaOFeAs single crystals cleave be- 
tween the La and As layers forming a polar surface. In 
contrast LiFeAs, of which clean large single- crystals are 
now available [8l4l0j. can be expected to cleave between 
adjacent Li layers and hence to produce a neutral sur- 
face. This is confirmed by results of density functional 
(DFT) calculations presented below. 

Although LiFeAs in bulk is found much more three 
dimensional (3D) compared to LaOFcAs in DFT calcu- 
lations [Ill4l3| . its surface electronic structure is much 
more bulk-like. Surface sensitive probes arc therefore ex- 
pected to directly yield information on the bulk behavior. 
Also, undoped stoichiometric LiFeAs is superconducting 
below 18 K and does not order magnetically. 13, 14 1 

Recent ARPES experiments on LiFeAs have revealed 
several interesting features. Although the low-energy 
electronic structure is qualitatively similar to other Fe- 
pnictides, the Fermi surface of LiFeAs is remarkably dif- 
ferent. It shows a strongly reduced tendency to (jr, 7r)- 
nestingflBl]. which is present in most of the new iron su- 
perconductors and often believed to be characteristic for 
all of them. As shown in Fig. [5j the radii of the hole and 
electron cylinders arc quite different here. Another ob- 
servation is the renormalization of the conduction bands 
by a factor up to 3 as compared to DFT calculations 
with an extended van Hove singularity (flat maximum of 
the h.2 band) in the center of the Brillouin zone. The 
four bands supporting the Fermi surface were found to 
be differently, though isotropically, gapped in the super- 



conducting state. [l5| These observations can have impor- 
tant implications for understanding the mechanism of su- 
perconductivity in the Fe-pnictidcs. Here we show that 
there is no surface driven electronic structure in LiFeAs 
and therefore surface sensitive measurements represent 
the bulk electronic structure as opposed to the situation 
in the 1111 (e.g. LaOFeAs) compounds. 

Non-spin-polarized scalar relativistic DFT total en- 
ergy and Kohn-Sham band structure calculations were 
performed with the full-potential local-orbital code ver- 
sion FPLO9.00 [lfij . All considered structures were fully 
relaxed with the generalized-gradient-approximation 



(GGA) functional [17] using the force tool of FPLO9.00. 
Forces were minimized to less then 5 • 10~ 3 eV/A for bulk 
and 10 _1 eV/A for slab calculations. This is necessary 
since no detailed surface structure is known from exper- 
iment. 

The calculations were carefully converged with the 
number of k-points, and finally 12,12,4 regular k-point 
meshes were used in the full Brillouin zone (BZ). The 
densities were converged down to 10~ 6 for the corre- 
sponding parameter in FPLO (rms(/9 out — p m )/u.c). Ta- 
ble U shows the results for the bulk lattice parameters 
compared to experiment [Hill 113- 
One can observe a decrease in volume of the cell, com- 
pared to experimental data, which is reflected in a c- 
lattice parameter reduction due to the softness of the 
lattice in z-direction (cf. LaOFeAs in Ref. [l]). This slight 
underestimation of the lattice parameters even in GGA 
calculations is a common feature for non-spin-polarized 
calculations in the pnictide family due to the strong cou- 
pling of magnetism to the lattice. However, the differ- 
ences remain acceptable although the deviation for the 
distance in z-direction between adjacent Fe and As layers 
dp c _As is approximately 3%. The local density approxi- 
mation (LDA), would result in even further reduced a- 
and c- values (results not shown), hence GGA was used 
throughout this work. 

After relaxation, the next step was to calculate the 
cleaving behavior of LiFeAs by the same approach as in 
Ref. [D. We had to use an orthorhombic Pmm2 space 
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a[A] 


c[A] 


ZAs 




dpc-As [A] 


experim. 


3.791 


6.364 


0.2635 


0.8459 


2.420 


GGA 


3.766 


6.158 


0.2729 


0.8325 


2.345 



Table I. Structure parameters of bulk LiFeAs; (1f c _as is the 
atom layer distance in ^-direction; experimental data taken 
from Ref . fil. 
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Figure 1. Inter layer distances of LiFeAs as a function of 
tensile strain in z-direction. The lines between data points 
are guides for the eye. 





bulk (4uc) 


3/1 slab 


3/3 slab 


bulk (6uc) 






Vac. 


Vac. 




Li/ As 


0.650 


0.496 


0.498 


0.650 


As/Fe 


1.398 


1.425 


1.424 


1.398 


Fe/As 


1.398 


1.396 


1.396 


1.398 


As/Li 


0.650 


0.634 


0.636 


0.650 


Li/Li 


2.062 


2.065 


2.066 


2.062 


Li/ As 


0.650 


0.643 


0.644 


0.650 


As/Fe 


1.398 


1.398 


1.399 


1.398 






Vac. 


Vac. 




Li 


-0.144 


-0.161 


-0.162 


-0.144 


As 


+0.296 


+0.426 


+0.426 


+0.296 


Fe 2 


-0.304 


-0.375 


-0.380 


-0.304 


As 


+0.296 


+0.269 


+0.267 


+0.296 


Li 


-0.144 


-0.155 


-0.152 


-0.144 


Li 


-0.144 


-0.146 


-0.144 


-0.144 


As 


+0.296 


+0.295 


+0.293 


+0.296 


Fe 2 


-0.304 


-0.304 


-0.297 


-0.304 



Table II. Upper part: Inter layer distances in A. Lower part: 
Electron excess per site. Left and right double column: see 
explanation in the text. In view of the mirror symmetry with 
respect to the central layer only half of a slab is represented. 



group instead of P4/nmm, in order to allow all atom 
layers to relax in z-direction, unconstrained by symmetry. 

In order to study the cleaving behavior we constructed 
a super-cell with two unit cells in z-direction, with lattice 
constant 2c. For c = Cbuik this just represents a bulk 
calculation. Now, we apply tensile strain in z-direction, 
by increasing c. We performed a set of calculations for 
c-values, which correspond to 0, 10, 20 and 30% strain. 

Fig. [T] shows the inter and intra unit-cell layer dis- 
tances of different atomic layers, such as Li-Li (straight 
and dashed lines), Fe-As (dotted line) and Li- As (dash- 
dot line) for this periodic superstructure of two unit cells 
in z-direction as a function of tensile strain. At large 
enough strain the original periodicity of period c in z- 
direction is broken by forming slits. By our adopted 
superstructure of two unit cells, slits form between ev- 
ery other Li-Li layer pair. (Inter indicates the dis- 
tance crossing the slit, when it is formed, while intra 
corresponds to distances within the two unit cell slab 
[Li\As\Fe 2 \As\Li\Li\As\Fe 2 \As\Li].) 

For small strain the whole slab is more or less uni- 
formly stretched with all the layer distances increasing 
compared to their bulk values. From the calculated forces 
of about 0.2 eV/A one infers a tensile stress in z-direction 
of 2.1 GPa for 10% tensile strain. For a strain between 
10 and 20%, LiFeAs cleaves between two Li layers in our 
calculations with a resulting neutral surface and conse- 
quently a non-polar slab, as we pointed out above (see 
also Ref. [H|). Cleavage here means that the layers are 
no longer uniformly stretched but that the two unit cells 
relax back to bulk like distances, while a slit (vacuum) 



opens up between two Li layers such that we have a re- 
peated slab setup with a slab thickness of two unit cells. 
In the doubled unit cell the interior Li-Li layer distance 
(termed intracell in Fig. [T]) relaxes back even below the 
bulk value. Also the topmost Li- As layer distance close 
to the slit reduces appreciably while the Fe-As layer dis- 
tance hardly changes. (It should be pointed out that our 
simulation of the cleaving process only provides its quali- 
tative behavior, because the chosen supercell may not be 
large enough to simulate true bulk like behavior in the 
interior layers.) 

To obtain a reliable lattice and electronic structure of 
a cleaved crystal surface it is necessary to choose a slab 
that is thick enough to pin the Fermi level of the slab to 
the bulk value and to choose a slit that is wide enough 
so that electronic wave functions from both sides of the 
slit do not overlap. (Due to the chosen mirror symmetry 
in the center of the slab (z-direction) there is no electric 
field across the slit.) It turns out that the Fermi level 
of even a slab thickness of only three unit cells hardly 
changes compared to the bulk, indicating that this thick- 
ness is already sufficient. To check convergence with the 
vacuum layer width, calculations were done with a slit 
width of about Cbuik (resulting in a periodicity 4c for the 
repeated three unit cells slab; termed 3/1 in Table IlT]) and 
with a slit width of about 3c (resulting in a periodicity 
of 6c for the whole repeated slab setup, termed 3/3). In 
order to have fc-integration errors under control, the bulk 
calculations (used for comparison in Table [TTJ) were done 
with supercells of the same periodicity in z-direction (4 
and 6 unit cells, respectively) allowing for the same k- 



Figure 2. Band structure of bulk LiFeAs, orbital projected 
onto the leading Fe, As and Li orbitals. 



mesh interpolation in both bulk and slab cases. The two 
different vacuum setups with their corresponding bulk 
calculations are represented by the left and right double 
columns in Table |TTJ 

When comparing the inner layer distances and layer 
charges of the slab with those of the bulk crystal (lower 
parts of both sub panels in Table [TTJ) , one sees that a 
slab of three unit cells thickness is already giving bulk 
representative interior layers. From the residual forces 
the numerical accuracy of the distances is estimated to 
be better than O.OlA, while the accuracy of the layer 
charges is better than 0.005 of the electron charge. Also 
the positions and charges between a calculation with a 
vacuum width of c and 3c barely differ, so that in what 
follows the 3/1 slab is used for the surface band structure 
analysis. 

As is seen from the upper parts of Table [HI there is a 
slight outward shift of the topmost As layer beneath the 
surface and an appreciable inward shift of the Li surface 
layer. Interestingly, this causes a charge redistribution 
between the As layer and the Fe layer below of about 
0.1 electron charge, but nearly no change of the Li layer 
charge. 

After having justified our model, we will inspect more 
closely Figs. [2] to HI which show results for the band- 
structure and Fermi surface (FS). In Fig. [2] we present 
the orbital projected band structure for the bulk calcu- 
lation. The orbital projection is obtained by projecting 
the Kohn-Sham wave functions onto the chemical orbitals 
of the FPLO basis. The line width of the colored lines 
in Fig. [2] indicates the orbital content/character in the 
wave function for every band. As is well known the Fe 
3d bands (shown in red) together with some hybridiza- 
tion of the As 4p states (green) dominate the entire low 
energy spectrum. This is responsible for the metallic be- 
havior of LiFeAs. 

However, it is important to note the strong z- 
dispersion of a Li band (shown in blue), coming from 
about 5 eV between M and T down to V (at 3c V) and 
going further down to -2 eV via several branches along 



Figure 3. Band plot of bulk LiFeAs for fc z = (black) and 
k z = 7r (red) and slab LiFeAs (green lines). 

T-M and T-Z-R-A, which couples to As at higher ener- 
gies and to Fe at lower energies. This effect is absent 
in all other families of iron pnictides or chalcogcnidcs 
i.e., 11 (e.tr FeSe), 1111 (e.g. LaOFeAs) or 122 (e.g. 
BaFe2As2)[ll(. In fact, this dispersive feature is respon- 
sible for the closed Fermi surface (FS) seen around T in 
Fig. 2] (left panel). Of course, in a slab model there is 
no fc z -dispersion and hence this FS is not closed in Fig. I 
(right panel). 

As seen in Fig. [31 all slab bands (thin green lines) in 
the vicinity of the Fermi level are within the fc z -dispersion 
range of the corresponding bulk bands (pairs of black and 
red lines). Hence no surface band is found in an energy 
window free of bulk bands. This is in contrast to the situ- 
ation found for LaOFeAs in RefH Furthermore, the bulk 
FS nesting is less pronounced than in LaOFeAs (larger 
difference between radii of electron and hole FSs), due 
to a larger k z -dispersion compared to LaOFeAs, which is 
mainly caused by the extended Li states. This could be 
related to the absence of magnetic order in this material. 

We note that Ref. [l9| calculated a smaller orbital 
weight ratio between As 4p and Fe 3d. Although both 
calculations differ somewhat in the Fe-As layer distance 
(here 2.345 A compared to 2.403 A of Ref. [H) this is 
mainly caused by different local orbital projections. 

Fig. Q] shows the calculated FSs of the bulk and the 
slab. As can be seen the fc^-dependence of the bulk FSs 
interpolates between the corresponding A^-dispersionless 
slab FSs, and no indication of surface specific features is 
found. Hence, surface sensitive experiments will repre- 
sent the bulk electronic structure. 

ARPES results for the FS are shown on Fig. [5] The FS 
radii are comparable with the calculated results at least 
in chosen directions while the shape of the FSs around T 
seems to be more complex in experiment. We infer that 
this finding, just as the band renormalization presented 
below, is intrinsic to the bulk electronic structure. 

In Table IIIII the calculated bulk Fermi radii (averaged 
over k z ) are compared to measured values (see Fig.[5|). 

In Fig. [H] we compare the experimental ARPES data 
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Figure 4. (color online) FS of LiFeAs. The Z-point is above 
r, the A-point above M and, the R-point above X. Left panel: 
bulk. Right panel: slab (no z-dispersion; 3D perspective for 
better visibility only). 
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Figure 5. (color online) FS of a LiFeAs from angle resolved 
photoemission. 



with the results of the calculations. The spectra were 
collected at BESSY. A detailed description of the exper- 
iment can be found elsewhere (Ref. [l5|). In the upper 
panel, we show the photoemission intensity and in the 
lower panel the corresponding results for the calculated 
bulk band structure. The photoemission data have been 
interpolated by a tight-binding (TB) fit (denoted by the 
dotted lines). Note, however, that there are uncertain- 
ties in TB fits, especially concerning band connectivities. 
The hole-like structures near the r-point are indicated 
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M 2 


Ti 


r 2 


experim. 
GGA 


0.15 
0.17 


0.18 
0.17 


0.24 
0.19 


w 0.07 
0.13 



Table III. Fermi radii in fc^-direction of different FS-pockets in 
units of — from experimental data and from the calculation. 
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Figure 6. (color online) Upper panel. Photoemission inten- 
sity recorded along the high symmetry directions in LiFeAs. 
Lower panel. Corresponding results for the calculated bulk 
bandstructure. 



by hi, li2 and I13, and electron- like structures near the 
corner of the BZ by ei and e2. These notations are also 
used in Fig. [S] (except for h 3 , see below). 

The energy scale in both panels is of the same width, 
which allows to estimate the band renormalization. The 
overall agreement is remarkable, taking into account that 
in both cases the low energy electronic structure is formed 
by five bands and the Fermi surface (FS) is of the same 
topology. Applying a renormalization factor to the cal- 
culations leads to even better, quantitative agreement. 
All characteristic dispersions and Fermi velocities are re- 
produced by the calculations. The experimental Fermi 
surface (Fig. [5]) consists of two hole-like FSs around the T- 
point and two electron- like FSs centered at the M-point . 
The strongly dispersive h 3 feature (seen in Fig. 16} barely 
contributes to the spectral weight at the Fermi level and 
hence is not indicated in Fig. [SJ The relative sizes of the 
Fermi surfaces differ slightly from the theoretical ones. 
Notably, the first two hole pockets at the zone center 
seem to be smaller in the ARPES data compared to the 
calculations. (The agreement can be improved by apply- 
ing a small shift of the chemical potential. This would 
also improve the agreement concerning the I13 feature at 
the r-point.) 

In summary, we have investigated the surface structure 
and surface effects on the electronic structure in LiFeAs. 
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The main finding is that due to the neutral cleaving of the 
compound between the Li-layers there are no significant 
surface effects on the electronic states. Consequently, 
LiFeAs is an ideal compound in the pnictide family for 
the study by ARPES and STM methods. 

The project was supported, in part, by the DFG under 
Grants No. KN393/4, BO 1912/2-1 as well as priority 
program SPP1458. 



[7] 



[8] 



[9] 
[10] 



* [hTeschrig@ifw-dresden.de; | http:/ /www.ifw-dresden.de/~/helmut| 
[1] H. Eschrig, A. Lankau, K. Koepernik, Phys. Rev. B 81, 
155447 (2010) 

[2] D. H. Lu, M. Yi, S.-K. Mo, A. S. Erickson, J. Analytis, 
J.-H. Chu, D. J. Singh, Z. Hussain, T. H. Geballe, R. R. 
Fisher and Z.-X. Shen, Nature 455, 81 (2008) 

[3] Ch. Liu, T. Kondo, A. D. Palczewski, G. D. Samolyuk, 
Y. Lee, M. E. Tillman, N. Ni, E. D. Mun, R. Gordon, 
A. F. Santander-Syro, S. L. Bud'ko, J. L. McChesney, 
E. Rotenberg, A. V. Fedorov, T. Valla, O. Copie, M. 

A. Tanatar, C. Martin, B. N. Harmon, P. C. Canfield, 
R. Prozorov, J. Schmalian and A. Kaminski, Physica C 
469, 491 (2009). 

[4] D. H. Lu, M. Yi, S.-K. Mo, J. G. Analytis, J.-H. Chu, A. 

S. Erickson, D. J. Singh, Z. Hussain, T. H. Geballe, I. R. 

Fisher and Z.-X. Shen, Physica C 469, 452 (2009). 
[5] Ch. Liu, Y. Lee, A. D. Palczewski, J.-Q. Yan, T. Kondo, 

B. N. Harmon, R. W. McCallum, T. A. Lograsso and A. 
Kaminski, larXiv:1006.0929V l 

[6] L. X. Yang, B. P. Xie, Y. Zhang, C. He, Q. Q. Ge, X. F. 
Wang, X. H. Chen, M. Arita, J. Jiang, K. Shimada, M. 
Taniguchi, I. Vobornik, G. Rossi, J. P. Hu, D. H. Lu, Z. 
X. Shen, Z. Y. Lu, D. L. Feng, larXiv:1006.1107i d 



[11] 
[12] 
[13] 
[14] 

[15] 

[16] 
[17] 
[18] 

[19] 



H. Liu, G. F. Chen, W. Zhang, L. Zhao, G. Liu, T.-L. 
Xia, X. Jia, D. Mu, Sh. Liu, Sh. He, Y. Peng, J. He, Zh. 
Chen, X. Dong, J. Zh ang, G. Wang, Y. Zhu, Z. Xu, Ch. 
Chen and X. J. Zhou. larXiv:0912.2838H (2010) 
Igor Morozov, Alexander Boltalin, Olga Volkova, Alexan- 
dre Vassiliev, Olga Kataeva, Ulrike Stockert, Mahmoud 
Abdel-Hafiez, Madeleine Fuchs, Hans- Joachim Grafe, 
Giinter Behr, Riidiger Klingeler, Sergey Borisenko, 
Sabine Wurmehl, Bernd Biichner, submitted to Cryst. 
Growth and Design 

Y. J. Song, J. S. Ghim, B. H. Min , Y. S. Kwon, M. H. 

Jung, J.-S. Rhvee. larXTv:1002. 22491 (2010) 

Joshua H. Tapp, Zhongjia Tang, Bing Lv, Kalyan Sas- 

mal, Bernd Lorenz, Paul C. W. Chu, and Arnold M. 

Guloy, Phys. Rev. B 78, 060505(R) (2008) 

H. Eschrig, K. Koepernik, Phys. Rev. B 80, 104503 

(2009) 

Takashi Miyake, Kazuma Nakamura, Ryotaro Arita, 
Masatoshi Imada, J. Phys. Soc. Jpn., 79, 044705 (2010) 
Hiroki Nakamura, Masahiko Machida, Tomio Koyama, 
Noriaki Hamada, J. Phys. Soc. Jpn. 78, 123712 (2009) 

C. W. Chu, F. Chen, M. Gooch, A.M. Guloy, B. Lorenz, 
B. Lv, K. Sasmal, Z.J. Tang, J.H. Tapp, Y.Y. Xue, Phys- 
ica C, 469, 326-331 (2009) 

S. V. Borisenko, V. B. Zabolotnyy, D. V. Evtushinsky, T. 
K. Kim, I. V. Morozov, A. N. Yaresko, A. A. Kordyuk, 
G. Behr, A. Vasiliev, R. Follath, B. Biichner Phys. Rev. 
Lett. 105, 067002 (2010) 

K. Koepernik and H. Eschrig, Phys. Rev. B 59, 1743 
(1999) 

J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

M. J. Pitcher, D. R. Parker, P. Adamson, S. J. C. Herkel- 
rath, A. T. Boothroyd, R. M. Ibberson, M. Brunelli and 
S. J. Clarke, Chem. Commun. 45, 5918 (2008) 

D. J. Singh, Phys. Rev. B 78, 094511 (2008) 



